Frustrated ferromagnetic spin-1/2 chain in a magnetic field: The phase diagram and 

thermodynamic properties 
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The frustrated ferromagnetic spin-1/2 Heisenberg chain is studied by means of a low-energy field 
theory as well as the density-matrix renormalization group and exact diagonalization methods. 
Firstly, we study the ground-state phase diagram in a magnetic field and find an 'even-odd' (EO) 
phase characterized by bound pairs of magnons in the region of two weakly coupled antiferromagnetic 
chains. A jump in the magnetization curves signals a first-order transition at the boundary of the 
EO phase, but otherwise the curves are smooth. Secondly, we discuss thermodynamic properties at 
zero field, where we confirm a double-peak structure in the specific heat for moderate frustrating 
next-nearest neighbor interactions. 



The physics of frustrated quantum spin systems is 
currently attracting large interest as exotic quantum 
phases may emerge. 1 Prominent examples are quantum- 
disordered ground states with different patterns of broken 
translational symmetry and quantum chiral phases (see, 
e.g., Ref.0). In addition, some frustrated systems have a 
large number of low- lying excitations, leading to unusual 
features in thermodynamic quantities. 

In one dimension, the paradigmatic model is the frus- 
trated spin-1/2 chain: 



H = VVis*; • + J2S1 ■ si 
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(1) 



si are spin s = 1/2 operators at site I, while h denotes a 
magnetic field. 

Much is known about the ground-state properties and 
the magnetic phase diagram of the frustrated antiferro- 
magnetic (AFM) chain with J 1; J 2 > 0. 2 We highlight 
the appearance of a plateau in the magnetization curve 
at magnetization M = 1/3 and the existence of an 'even- 
odd' (EO) region at small J\ with spins flipping in pairs 
in a magnetic field h£ 

Relatively little attention has been paid to frustrated 
ferromagnetic (FM) chains, i.e., J\ < and J2 > 0, until 
the recent discovery of materials described by parame- 
ters with this combination of signs. We mention in par- 
ticular Rb2Cu2Mo 3 0i2 which is believed to be described 
by Ji ps — 3 and LiCuV04 which lies in a differ- 
ent parameter regime with J\ sa —0.3 J2? In both cases, 
the saturation field h aa t is within experimental reach. 
A recent transfer-matrix renormalization group (TMRG) 
studj£ of the thermodynamics of Eq. (JTJ was motivated 
by the experimental results for Rb2Cu2Mo30i2- 

In this paper we study the zero-temperature phase di- 
agram in a magnetic field and the thermodynamics of 
Eq. {TJ at zero field. The former is obtained by a combi- 
nation of a low-energy field theory and the density ma- 
trix renormalization group (DMRG) method^ while the 



latter is computed by exact diagonalization (ED). We 
develop a minimal effective field-theory description for 
the region of small J± and h < h sa t and predict the ex- 
istence of an EO phase. Note that at J\ — — 4 J2, the 
system undergoes a transition to a FM ground stated 
The field theory predictions are verified by our DMRG 
results. Further, our ground-state phase diagram differs 
qualitatively from recent mean- field predictions A In our 
study of thermodynamic properties^ we focus on the 
example of J± = —3 J2 and present data for system sizes 
up to N = 24 sites. The specific heat of LiCuV04 will 
be discussed elsewhere iii 

First we discuss an effective field theory describing the 
long wavelength fluctuations of Eq. Q in the limit of 
strong next-nearest neighbor interactions J2 ^> | J\ \ . 

Just below the saturation field, the problem can be 
mapped onto a dilute gas of bosons^ This mapping, 
which is asymptotically exact in the subspace of two 
magnons, shows that magnons bind in pairs for any 
J\ < 0. Although the two-magnon state is not always re- 
alized as a ground state in a magnetic fieldji^ Chubukov 12 
found that in this subspace and for — 0.38Ji < J2, 
the ground-state momentum is commensurate while for 
— 0.25Ji < J2 < — 0.38Ji, it becomes incommensu- 
rate. Based on the discontinuous nature of the change 
of momentum for the lowest two-magnon bound state, 
Chubukov further predicted a first-order phase transition 
between a chiral and a dimerized nematic-like phase. 

Apart from the issue of the two-magnon states being 
realized as ground states, the mapping onto a dilute gas 
of bosons is controlled just near the saturation field h « 
hs&t ■ We apply a complementary bosonization procedure 
which is controlled for h < h sat and confirm that the 
hallmark property of the commensurate region - pair- 
binding of magnons - is universal and extends well below 
the saturation field. A good starting point is the limit 
of J2 ^> |Ji| and a finite magnetic field h 7^ 0~ In 
this limit, the system may be viewed as two AFM chains 
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subject to an external magnetic field and weakly coupled 
by the FM zig-zag interaction J\. It is well known that 
the low-energy effective field theory for a single isolated 
spin-1/2 chain {J\ = 0) in a uniform magnetic field is the 
Tomonaga-Luttinger liquid: 14 



v 

H =2 



dx{^(d x 0) 2 + K(d x 6) 2 }. 



(2) 



Above we have introduced a compactified scalar bosonic 
field <f> and its dual counterpart 8, with [<j)(x), 9(y)] = 
iQ(y — x), where Q(x) is the Heaviside function. 

The Luttinger liquid (LL) parameter K(h) and spin- 
wave velocity v(h) can be related to microscopic parame- 
ters of the lattice model J2 and h using the Bethe-ansatz 
solution of the Heisenberg chain in a magnetic fieldi 1 ^* 1 ^ 
We recall here that K(h) increases monotonically with 
the magnetic field from K(h = 0) = 1/2 to the uni- 
versal free-fermion value K = 1 for h approaching the 
saturation field h sa t — 2J2. The Fermi wave vector 
kp — -| (1 — M) is determined by the magnetization M. 
Note that we normalize the magnetization to M = 1 at 
saturation, i.e., M = 2S Z /N with S z = J2i s i- 

Now we perturbatively add the interchain coupling 
term to two chains, each of which described by an ef- 
fective Hamiltonian of the form Eq. J2J and fields </>j, 
i = 1, 2. For convenience, we transform to the symmetric 
and antisymmetric combinations <j>± = (4>i i^)/v* arL d 
9± = (9i ± 62) /V%- In this basis and apart from terms 
H.^ of the form © , the effective Hamiltonian describing 
low-energy properties of Eq. Q contains a single relevant 
coupling with the coupling g\ oc J\ -C v: 

H e g = H$ + Hq + pi J dx cos {kp + V8n(j)^) , (3) 

and the renormalized LL parameter: 



~i 1 1 1 1 1 r 




Further 

phases? 
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K- =K{h) {1 + J x K{h)/[irv(h)\} . (4) 

The Hamiltonian Q yields the minimal effective low- 
energy field theory describing the region J2 3> | J\ \ of the 
frustrated FM spin-1/2 chain for M^O. The relevant 
interaction term cos opens a gap in the <f>- sector. 

Since sf +1 — sf ~ d x 4>_, relative fluctuations of the two 
chains are locked, leading to pair-binding of magnons. 
These bound pairs of magnons themselves are gapless, 
since sf +1 + sf ~ d x 4>+- This phase was observed for an 
AFM Ji in Ref.0 and dubbed the 'EO phase'. 

In addition, we confirm this picture numerically. The 
simplest possible lattice model that is described by a low- 
energy effective Hamiltonian of type © is a spin ladder 
with a dominant biquadratic leg-leg interaction^ We 
compute the magnetization curve with DMRG (results 
not presented here) and verify that only even magnon 
sectors are realized as ground states for all fields h > 0. 

Equation Q shows that the LL parameter K_ de- 
creases with an increasing absolute value |Ji| of the FM 
interchain coupling, in contrast to an AFM coupling. 



FIG. 1: (Color online) (a), main panel (inset): Magnetization 
curve M(h) for Ji = — J 2 (Ji = — 2.5J2). The horizontal 
dotted line marks M = 1/3. (b): M(h) for Ji = -3 J 2 . (c): 
Magnetic phase diagram of the frustrated FM chain. The dot- 
ted line (with stars) marks the first-order transition between 
the EO phase and the AS Z = 1 region, while the line h = ht 
(dashed, triangles) separates the AS 12 = 1 region from the 
AS Z > 3 part. Uncertainties of the transition lines, e.g. due 
to finite-size effects, should not exceed the size of the sym- 
bols. The fields /ij ump and ht were extracted from iV = 156 
sites (stars) and JV = 120 sites (triangles], respectively. The 
dashed, vertical line is the result of Ref. Ha (J2 « 0.38Ji). 

However, the bosonization procedure becomes inappli- 
cable once K_ vanishes. This signals an instability of 
the EO phase when increasing \J\\. Moreover, we con- 
clude that for FM J\, the EO phase extends up to the 
saturation field h sat (since K- < 1 for J\ < such that 
there is always a relevant coupling in the antisymmetric 
sector), in contrast to the AFM case. 

To check this scenario and to determine the phase 
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boundaries, we perform DMRG calculations for up to 
156 sites imposing open boundary conditions. The finite- 
system algorithm! is used and we keep up to 350 states. 
DMRG gives direct access to the ground-state energies 
Eo(S z ,h = 0) at zero magnetic field in sub-spaces la- 
beled by S z . After shifting the ground-state energies 
E (S z ,h = 0) by a Zeeman term through E (S z ,h) — 
Eo(S z ,h — 0) — hS z , it is straightforward to construct 
the magnetization curve. 

We start the discussion from the limit |Ji| < J2. 
The magnetization curves for J\ — —Ji are shown in 
Fig. H ( a )- I n particular, we verify the pair- binding of 
magnons predicted above: in a wide parameter range in 
the magnetic phase diagram (h vs J\), the magnetization 
changes in steps of AS Z = 2. This can be observed even 
on systems as small as N — 24, while for an AFM inter- 
chain coupling Ji > 0, the formation of bound states was 
only reported for long chains^ 

From the inset of Fig. ^ ( a )> which shows data for 
J 1 = — 2.5 J 2, we conclude that a second phase emerges 
at lower fields, signaled by a change of the magnetization 
steps from AS* 2 = 1 to 2 at h = ft.j U mp- This transition is 
first order. 

In contrast to the frustrated antiferromagnetic chain^ 
no indications of a M = 1/3 plateau are found. We find 
that the width of the 1/3 plateau as seen on finite systems 
scales to zero with 1/N. 

The magnetization curve exhibits further features 
when Ji approaches the transition to the FM regime, 
occurring at J\ = — 4 J2. The main observations from 
Fig-D(b), which shows M(h) for J\ = — 3 J2, and ad- 
ditional data not displayed in the figures, are as follows. 
Below saturation, the steps in M(h) are of size AS Z — 3 
[see, e.g., the case of J\ = —3 J2 in Fig.Q](b)]. Upon de- 
creasing Ji — * — 4 J2, the magnetization curve becomes 
very steep below saturation and the steps of AS Z may 
even be larger than 3. For instance, we find steps of 
AS* 2 = 4 for Ji = —3.75 and N — 60 below saturation. 
Nevertheless, the asymptotic behavior of M(h) close to 
the saturation field for Ji = —3 J2 [see Fig.Q](b)] is con- 
sistent with a standard square-root singularity in M(h) 
(see Ref . IT3I and references therein). 

Our main findings for the magnetic phase diagram are 
displayed in Fig.n(c), where h is normalized by /i sat >i2ii£ 
The largest part of the phase diagram belongs to the EO 
phase, while the transition to the region with AS" 2 = 1 is 
first order. The position of this line, i.e., h;„m r,, (d otted, 
with stars) is consistent with results of Ref. [12 in the 
high- field limit, but the transition takes place at lower 
h/h S nt for smaller J\. For larger J\ and fields h > h t , 
a third region emerges, characterized by AS Z ~ 3. Just 
as for AFM Jir^ one may speculate about chiral order 
in some of these regions as well as additional phases, but 
substantially larger system sizes might be needed to fully 
reveal the nature of this part of the phase diagram. 

Next we discuss thermodynamic properties concentrat- 
ing on h = 0. We perform full diagonalizations to ob- 
tain all eigenvalues and then use spectral representations 
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FIG. 2: (Color online) Magnetic susceptibility (top panel), 
specific heat (middle panel) and entropy per site (bottom 
panel) for N = 16, 20, and 24 at Ji = -3J 2 , h = in 
comparison to a FM chain with J2 = and iV = 20. Middle 
panel, inset: specific heat at low temperatures for Ji = —3 J2. 



to compute thermodynamic quantities, as described in 
some detail for the entropy in Ref. In order to ren- 
der the Hamiltonian Q translationally invariant, we now 
impose periodic boundary conditions. After symmetry 
reduction, the biggest matrices to be diagonalized for 
— 24 are of complex dimension 81 752. In such high 
dimensions, we use a custom shared memory parallelized 
Householder algorithm, while standard library routines 
are used in lower dimensions. 

Figure [5] shows results at J2 = — Ji/3 and h — 
for rings with N = 16, 20, and 24. This ratio 
of exchange constants is close to values suggested for 
Rb 2 Cu2Mo30i2)^ and the phase diagram in a magnetic 
field promises interesting properties in this parameter 
regime. Both the magnetic susceptibility x an d the 
specific heat C have a maximum at low temperatures, 
namely for N — 24 at T w 0.04 |Ji| in the case of \ 
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and T w 0.023 |Ji| in the case of C. While these low- 
temperature maxima are affected by finite-size effects, 
the dependence on N is negligible at higher tempera- 
tures. The specific heat exhibits a second broad maxi- 
mum around T rs 0.67 |Ji|. Such a double-peak struc- 
ture in the specific heat has already been observed for 
J2 = —0.3 Ji on a finite lattice with N = 16 sites^Si and 
by TMRG at J 2 = -0.4 J v 6 Note that the results for C 
of Ref. are restricted to temperatures T > 0.013 |Ji| in 
this parameter regime, and the TMRG method might be 
plagued by convergence problems at low temperatures. 
Despite the finite-size effects in our data at low tempera- 
tures we can clearly resolve the low-temperature peak in 
C (see inset of the middle panel of Fig. |5J). 

Our results for \ (top panel of Fig. |2J differ qualita- 
tively from those obtained for J2 = — 0.3 Ji and N = 16 
in Ref. l2llbv ED. In particular, we find a singlet ground 
state for all periodic systems with | Ji | < 4 J2 investi- 
gated, in contrast to Ref. |^ However, we do find good 
agreement with the more recent TMRG results for x~ 

It is not entirely trivial to separate the low- and high- 
temperature features in C and \ m to FM and AFM ones. 
Let us compare the case J 2 = — </i/3 with an unfrustrated 
FM chain (Fig. [5] includes results for J 2 = 0, J\ < and 
N = 20). In both cases, there is a broad maximum in 
C at high temperatures, although numerical values are 
different. Concerning the low-temperature peaks in \ 
and C, note that, for J x = -3J 2 , the FM s = N/2 
multiplet is located at an energy of about N | J\ | /40 above 
the s — ground state. Since this energy scale roughly 
agrees with the temperature scale of the low-temperature 
maxima, it is conceivable that they correspond to FM 
fluctuations above an AFM ground state. 

Finally, we note that the entropy of the frustrated FM 
chain (J2 = — Ji/3) is larger than that of the simple 
FM chain (J2 = 0) over a wide temperature range (see 
bottom panel of Fig. |2J . Only for very low temperatures 



the FM ground state leads to a bigger entropy for J 2 = 0. 

To summarize, we have studied the ground state phase 
diagram of a frustrated FM chain in a magnetic field 
and found an EO phase characterized by bound pairs of 
magnons. The boundary of this phase appears to be first- 
order and terminates for h — ► h saX at J2 ~ — 0.38 
At larger FM | Ji | , changes in the step height of the mag- 
netization curves signal the presence of further phases, 
which need to be studied in more detail. It would also be 
desirable to better understand the low-lying excitations 
in the different phases and to compare to the case of the 
frustrated antiferromagnetic chain^ Our phase diagram 
differs substantially from recent mean-field predictions^ 
In particular, our DMRG data exhibit a smooth tran- 
sition to saturation for any J\ > — 4 J2, in contrast to 
previous studies This observation may also be rel- 
evant for the transition to saturation in the frustrated 
square lattice ferromagnet, 22 The parameters relevant to 
LiCuVO^ lie well inside the EO phase where the theo- 
retical magnetization curves are completely smooth. 

Furthermore, we have discussed thermodynamic 
properties^ The most prominent feature for J 2 — 
— Ji/3, h = is a double-peak structure in the spe- 
cific heatASI The excitation spectrum is not reflected 
directly in thermodynamic quantities, but microscopic 
probes such as neutron scattering or nuclear magnetic 
resonance should be able to differentiate between gapped 
AS Z = 1 excitations and gapless AS Z — 2 excitations. 
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